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Abstract 

Background and Purpose: Currently there are more and more studies on the association between short-term effects of 
exposure to particulate matter (PM) and the morbidity of stroke attack, but few have focused on stroke subtypes. The 
objective of this study is to assess the relationship between PM and stroke subtypes attack, which is uncertain now. 

Methods: Meta-analyses, meta-regression and subgroup analyses were conducted to investigate the association between 
short-term effects of exposure to PM and the morbidity of different stroke subtypes from a number of epidemiologic studies 
(from 1997 to 2012). 

Results: Nineteen articles were identified. Odds ratio (OR) of stroke attack associated with particular matter ("thoracic 
particles" [PM qo ]<10 urn in aerodynamic diameter, "fine particles" [PM 25 ]<2.5 |xm in aerodynamic diameter) increment of 
10 ug/m 3 was as effect size. PM 10 exposure was related to an increase in risk of stroke attack (OR per 10 ug/m 3 = 1.004, 
95%CI: 1.001-1.008) and PM 25 exposure was not significantly associated with stroke attack (OR per 10 ug/m 3 = 0.999, 
95%CI: 0.994-1.003). But when focused on stroke subtypes, PM 2 5 (OR per 10 ug/m 3 = 1.025; 95%CI, 1.001-1.049) and PM 10 
(OR per 10 pg/m 3 = 1.013; 95%CI, 1.001—1.025) exposure were statistically significantly associated with an increased risk of 
ischemic stroke attack, while PM 25 (all the studies showed no significant association) and PM 10 (OR per 10 ug/m 3 = 1.007; 
95%CI, 0.992—1.022) exposure were not associated with an increased risk of hemorrhagic stroke attack. Meta-regression 
found study design and area were two effective covariates. 

Conclusion: PM 25 and PM 10 had different effects on different stroke subtypes. In the future, it's worthwhile to study the 
effects of PM to ischemic stroke and hemorrhagic stroke, respectively. 
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Introduction 

Many studies regarded air pollution exposure as an important 
factor of hospitalization and mortality worldwide. PM, playing an 
important role in pollutants of major public health concern, had 
been confirmed that it could impair the respiratory and 
cardiovascular system through a series of changes in autonomic 
nervous system activity [1] and systemic inflammation [2], giving 
rise to alterations in oxidative stress [3,4], hematologic activation 
[5] and vascular endothelial dysfunction [6]. Most researches 
regarded PM 10 and PM 2 5 as major harmful PMs. 

Nevertheless, short-term effects of PM exposure on cerebral 
vessels were uncertain. Wordley et al. [7] and Tsai et al. [8] found 
that PM 10 was associated with daily stroke attack positively. While, 
in the works of Chan et al. [9], Henrotin et al. [10] and Andersen 
et al. [11], no significant association was demonstrated between 



PMio and hemorrhagic stroke attack. Similarly, analyses on the 
relationship between PM 2 5 and stroke attack also appeared to 
divergent results. Villeneuve et al. [12] found PM 2 5 exposure 
wasn't related to an increased risk of ischemic stroke attack (OR 
per 10 u.g/m 3 = 1.052, 95% CI: 0.996-1.160), while Wellenius et 
al. [13] found a positive association between PM 2 5 exposure and 
ischemic stroke attack (OR per 10 u.g/m 3 = 1.278, 95%CI: 
1.079-1.525). 

Our previous research focused on the association between PM 
exposure and stroke attack in two study designs (time-series design 
and case-crossover design), and the result indicated that the effects 
of PM to stroke attack varied in different study designs [14]. 
However, in addition to study design, there were still many other 
covariates (e.g. age, gender, economic condition, area, lags times, 
historical disease and temperature) among studies, which could 
influence the results. Of special interest was that whether PM can 
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act differently on different stroke subtypes. So in this article we 
determined to do meta-analyses, meta-regression and subgroup 
analyses of association between PM and different stroke attack. 

Methods 

1. Studies selection 

We identified studies published in English and Chinese up to 
March 2013, by literature search using PubMed, Web of Science, 
MEDLINE, Google Scholar, China National Knowledge Infra- 
structure (CNKI) and reference lists of relevant articles. Search 
terms included "Air Pollution/Particulate Matter" plus "Cardio- 
vascular disease/Stroke", besides, key terms "hospitalization/ 
Hospital Administration/Emergencies/Morbidity", "cardiovascu- 
lar disease" were used to enlarge the searching range. We chose 
ICD9: 430-438 or ICD10: 160-169 as the definition of "stroke" or 
"cerebrovascular disorders", ICD9: 430-432 or ICD10: 160-162 
for "hemorrhagic stroke", ICD9: 433-434 or ICD10: 163-166 for 
"ischemic stroke". 

Eligible studies were selected by two reviewers (X.L., J.S.) 
independendy according to following inclusion criteria: (1) The 
outcome focused on the effect of PM to stroke or cerebrovascular 
disease (2) Published full-text articles (3) focused on PM 10 and/or 
PM 2 .5 (4) Studies with similar effects [e.g. risk ratios (RR), 95% 
CIs] that could approximate ORs. The exclusion criteria were: (1) 
Duplications (2) Reviews or Meta-analysis (3) Long-term effects 
articles (4) air pollution from industrial and occupational 
environment (5) Articles that did not provide calculable or 
reported ORs and 95% CIs. The two reviewer reached consensus 
on the eligibility of each study. When there was disagreement, a 
third reviewer (G.C.) took the final decision. 

Using a standardized form, data from eligible studies were 
extracted by two reviewers (X.L., J.S.) independently, From each 
study, we collected name of first author, year of publication, 
number of participants, PM, country, age, gender, types of stroke, 
lags of air pollutants' concentrations, study design and ORs or 
RRs with 95% CI. In Meta-regression, covariates were gender, 
area, lag times, study design and research period. Gender was 
represented as female, male and the whole population. Areas 
included Asia, Europe, North America and Oceania. Lag times 
were the same day on stroke attack and the previous 1, 2, 3, 4, 5 
day. Study designs were classified into time-series and case- 
crossover designs. Research period was divided into early stage 
(1992-1997), middle stage (1997-2004) and late stage (2004- 
2009). Given there were 5 covariates (gender, area, lag times of air 
pollutants' concentrations, study design and research period), we 
used Meta-regression model to detect any possible influence 
factors (P<0.05). 

2. Statistical analysis 

Pooled ORs of PM with 95% CI for stroke attack were 
calculated by using the fixed or random effects Meta-analysis of 
with 0_ and I 2 statistics given as the chosen measure of 
heterogeneity (The null hypothesis of this test is homogeneity). 
The 0_ and I 2 statistics were used to assess heterogeneity, where 
P<0.05 or I 2 >50% were considered as significant heterogeneity. 
We presented random effects pooled estimates when heterogeneity 
was detected; otherwise, we used the fixed model. We also 
produced forest plots to show ORs from each of the individual 
studies included in the meta-analyses and the estimation of the 
pooled OR. The sizes of the markers of each OR in the plots 
represent the relative weight each study contributed to the pooled 
estimation. We assessed publication bias visually through funnel 
plots and a weighted Egger's test. We also performed sensitivity 



analyses, whereby each article was omitted in turn, recalculating 
the pooled estimates under extreme conditions. Moreover, Meta- 
regression analysis was performed to figure out whether the 
association between PM and stroke attack was influenced by 
covariates. With a positive Meta-regression coefficient presented 
(P<0.05), we could recognize the influence of the given factor. All 
analyses were performed using software STATA version 12.0 
(StataCorp LP, College Station, TX, USA). 

Results 

1. Study selection and data extraction 

A total of 107 potentially relevant researches were identified by 
searching electronic databases and reference lists. 19 full-text 
articles were eligible for inclusion in this analysis and data were 
extracted from these studies [7-13,15-26]. The details of the 
selection process were presented in a flow chart in Figure 1. 
Publication years ranged from 1997 and 2012. In total, 9 countries 
were involved including Australia, Canada, China, Denmark, 
Finland, France, Italy, UK and US. 14 articles focused on PM 10 
and 9 were about PM 2 5 . In our previous research [14], we only 
extracted adjusted maximum effective value in each study. 
However, in this article, we took advantage of all the effect values 
(33 for PM 2 5 , 68 for PM 10 , 29 for hemorrhagic stroke, 42 for 
ischemic stroke) that fulfilled our study aim. The basic overview of 
the 19 included articles is given in Table 1. 

2. Meta-analysis of different PMs exposure to different 
stroke types 

2.1 Effects of PMj^ PM ]0 exposure to stroke 
attack. There were nine articles (containing 33 studies) that 
referred to the association of PM 2 5 and stroke attack. Since 
heterogeneity existed among studies (Q,= 67.09, P<0.05), the 
random effect model was conducted to calculate a pooled OR with 
95%CI. The Meta-analysis results indicated that PM 2 3 exposure 
wasn't related to an increased risk of stroke attack (OR per 10 ug/ 
m 3 = 0.999, 95%CI: 0.994-1.003), Forest plot and Funnel plot 
were shown in Figure 2 and Figure 3(A), respectively. Egger's test 
(t—0.98, P— 0.336) didn't find evidence of publication bias. 
Sensitivity analysis suggested that no individual study significandy 
affected the pooled effect size, indicating that the results for PM 2 5 
and daily stroke attack were statistically robust (Table SI in File 
SI). 

Fourteen articles (containing 68 studies) were included. The 
heterogeneity was significant (Q_= 223.25, P<0.05). With the 
random effect model, we pooled all 68 studies into the meta- 
analysis and found PMjo exposure was statistically significandy 
associated with an increased risk of stroke attack (OR per 10 ug/ 
m 3 = 1.004; 95%CI, 1.001-1.008). Forest plot was shown in 
Figure 4. Funnel plot was shown in Figure 3(B). Egger's test 
supported that publication bias was unlikely (f=0.80, P= 0.427). 
Sensitivity analysis showed that results for PM 10 and stroke attack 
were not robust to the inclusion of Tsai study [8] and Vidale study 
[25] (Table SI in File SI). 

2.2 Effects of PM 2 .5^ PM 10 exposure to hemorrhagic stroke 
attack. Since the effect values of all studies included were not 
statistically significant, so it had no sense to do the Meta-analysis. 
As a result, we could not just simply determine the association 
between PM 2 5 exposure and hemorrhagic stroke attack. 

We conducted meta-analysis for nineteen combinative effects of 
PM 10 and hemorrhagic stroke attack. Heterogeneity was detected 
(Q,= 39.82, P= 0.002) through heterogeneity test. With random 
effect model, we found no evidence for the association between 
PM 10 exposure and hemorrhagic stroke attack (OR per 10 ug/ 
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Figure 1. Flow chart of the selection process. 

doi:1 0.1 371 /journal.pone.0095682.g001 
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Figure 2. Forest plot of ORs for the association between PM 2 . 5 and stroke attack. OR indicates odds ratio; CI, confidence interval. 
doi:1 0.1 371 /journal. pone.0095682.g002 



m 3 = 1.007; 95%CI, 0.992-1.022). Forest plot was shown in 
Figure 5, Funnel plot was shown in Figure 3(C). Egger's test 
showed that publication bias was unlikely in the meta-analysis on 
association between PM 10 and hemorrhagic stroke attack (t= 0.71, 
P— 0.487). Sensitivity analysis suggested the results for PM 10 and 
hemorrhagic stroke attack were statistically robust (Table SI in 
File SI). 

2.3 Effects of PM 2i5l PM 10 exposure to ischemic stroke 
attack. Heterogeneity was observed among five articles 
(Q = 24.00, P= 0.031). With random effect model, PM 25 exposure 



was related to the risk of ischemic stroke attack (OR per 1 0 ug/ 
m' 5 = 1.025; 95%CI, 1.001-1.049). Forest plot was shown in 
Figure 6. Funnel plot was shown in Figure 3(D). Egger's test 
(f=3.71) P— 0.003) showed publication bias existed among 
studies. Sensitivity analysis showed that results for PM 2 .5 and 
ischemic stroke attack were not robust to the inclusion of Lisabeth 
study [22], Villeneuve study [12] and Wellenius study [13] (Table 
SI in File SI). 

Significant heterogeneity existed among six articles included 
(Q_= 98.04, P= 0.000). Association was demonstrated between 
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Funnel plot with pseudo 95% confidence limits 
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Figure 3. Funnel plots of Meta-analysis in different particular matters to different stroke types. OR indicates odds ratio. (A) Funnel plot 
of Meta-analysis of PM 2 .s exposure to stroke attack. (B) Funnel plot of Meta-analysis of PM, 0 exposure to stroke attack. (C) Funnel plot of Meta-analysis 
of PM, 0 exposure to hemorrhagic stroke attack. (D) Funnel plot of Meta-analysis of PM 2 .5 exposure to ischemic stroke attack. (E) Funnel plot of Meta- 
analysis of PM q0 exposure to ischemic stroke attack. 
doi:1 0.1 371 /journal.pone.0095682.g003 



PMio exposure and ischemic stroke attack (OR per 10 |ig/ 
m 3 = 1.013; 95%CI, 1.001-1.025) using random effect model. 
Forest plot was shown in Figure 7. Funnel plot was shown in 
Figure 3(E). Egger's test = — 1.61, P= 0.120) indicated no 
publication bias existed among studies of association between 



PMio and ischemic stroke attack. Sensitivity analysis showed that 
results for PM 10 and ischemic stroke attack were not robust to the 
inclusion of Andersen study [1 1] and Vidale study [25] (Table SI 
in File SI). 
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Figure 4. Forest plot of ORs for the association between PM 10 and stroke attack. OR indicates odds ratio; CI, confidence interval. 
doi:1 0.1 371 /journal.pone.0095682.g004 
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Figure 5. Forest plot of ORs for the association between PM 10 to hemorrhagic stroke attack. OR indicates odds ratio; CI, confidence 
interval. 
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3. Results of Meta-regression of different PMs exposure 
to different stroke types 

We had detected heterogeneity among studies about PM 2 5 
exposure to stroke attack, PM 10 exposure to stroke attack, PM 10 
exposure to hemorrhagic stroke attack, PM 25 exposure to 
ischemic stroke and PM 10 exposure to ischemic stroke. Meta- 
regressions were done in these studies to detect the influence 
factors. (Table S2 in File SI). 

When studying PM 2 .5 exposure to stroke attack, we found study 
design (coefficient 0.032, 95% CI 0.002 to 0.062, P= 0.035) the 
influence factor. Among studies of PM 10 exposure to stroke attack, 
area (coefficient 0.007, 95% CI 0.002 to 0.012, P= 0.007) were 
found to be the possible influence factor. (Table S2 in File SI). 

4. Subgroup analyses 

Results above supported that design was the main covariate of 
studies about PM 2 5 exposure to stroke attack, so we divided these 
studies into subgroups (time-series group and case-crossover 
group) according to study design. Case-crossover group revealed 
the positive association between PM 2 5 exposure and stroke attack 
(OR per 10 ug/m 3 = 1 .029; 95% CI, 1 .003- 1 .055). And area was 
found to be the main covariate of studies about PM 10 exposure to 



stroke attack, subgroups (Asia, Europe, North America and 
Oceania) analysis showed that in Asia and Europe, PMjo was 
associated with stroke attack (Table 2). 

Discussion 

The pooled results of meta-analyses showed that PM 10 exposure 
was related to an increase in risk of stroke attack (OR per 1 0 ugV 
m 3 = 1.004, 95%CI: 1.001-1.008) and PM 25 exposure was not 
significandy associated with stroke attack (OR per 10 pig/ 
m 3 = 0.999, 95%CI: 0.994-1.003). So it's meaningful to explore 
the effects of PM 10 and PM 2 5 to stroke attack, separately. 

However, as we all know, ischemic stroke and hemorrhagic 
stroke are two different types of stroke. The mechanisms differ 
from each other. So the PM may affect different strokes in 
different ways. Thus, it is valuable to explore the association 
between PM and ischemic stroke or hemorrhagic stroke, 
respectively. This was the first Meta-analysis that focused on the 
association between particular matter and different stroke subtypes 
attack. Our research showed that PM 2 5 (OR per 10 ugV 
m 3 = 1.025; 95%CI, 1.001-1.049) and PM 10 (OR per 10 u.g/ 
m 3 = 1.013; 95% CI, 1.001-1.025) exposure were statistically 
significandy associated with an increased risk of ischemic stroke 
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Figure 6. Forest plot of ORs for the association between PM 2 . 5 and ischemic stroke attack. OR indicates odds ratio; CI, confidence interval. 
doi:1 0.1 371 /journal.pone.0095682.g006 



attack, while PM 2 .5 (all the studies showed no significant 
association) and PM ln (OR per 10 ug/m 3 = 1.007; 95%CI, 
0.992~ 1.022) exposure were not associated with an increased 
risk of hemorrhagic stroke attack. W e inferred that this phenom- 
enon might be caused by the different formation mechanisms of 
different stroke types. For the biological mechanisms of PM to 
cerebrovascular disease are not clearly clarified at present, we can 
just supposed that PM's effect to ischemic stroke attack might be 
regulated by systemic inflammation and the activation of 
coagulation system, which leading to atherosclerosis, vasoconstric- 
tion, increase of fibrinogen and acceleration the formation of acute 
thrombus. And PM's effect to hemorrhagic stroke attack might be 
caused by systemic inflammation, giving rise to alterations in 
oxidative stress, vascular endothelial damage and rupture of 
plaque. According to our results, it might be easier for PM to 
activate coagulation system and constriction of blood vessels than 
simply destroy the vascular endothelial. But it was just a hypothesis 
that needed to be tested. 

Heterogeneity in the studies we retrieved could come from 
inherent differences between study settings, as well as from 
differences in age, gender, area, lag times of air pollutants' 
concentrations, study design and so on. With respect to these 
influence factors, we used randomized effect model to minimize 
the heterogeneity and do Meta-regression to find the possible 
covariates. As a result, we only found design and area were the 
influence factors in the studies. 



Sensitivity analysis revealed that results for PM and ischemic 
stroke were not robust, but results for PM and hemorrhagic were 
more robust than the former. Except for the studies of PM2.5 and 
ischemic stroke, all the studies showed a symmetric inverted funnel 
shape, which indicated publication bias was unlikely. Besides, 
Egger's test, of which the results support funnel shape, was used to 
quantitatively assess the publication bias. 

In this research, the following issues needed to be considered: (1) 
there were limitations for selecting studies. First, Meta-analysis 
articles [27-30] were not included in our research, because the 
methods authors used were different from ours, and they only 
provided pooled effects of PM exposure to stroke attack. Second, 
long-term effects articles [31-33] were not included, for there 
existed more covariates [e.g. body-mass index (BMI), smoking 
status, blood pressure, educational level, household income and 
historical disease] in long-term effect articles. Third, articles [34— 
37] proving quantitative relationship between PM and stroke 
attack but didn't provide calculable or reported ORs and 95% CIs 
were removed, too. Hence, removing these three kinds of articles 
might decrease bias but could lose some evidence of the 
association between PM and stroke attack. (2) Measurement of 
air pollution exposure varies within and between studies. The air 
pollution monitor apparatus itself had measurement error and was 
different from study to study. Besides, it was certainly known that 
personal pollution exposure levels were very different from those 
measured at a nearby fixed monitoring station [38]. Thus 
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Figure 7. Forest plot of ORs for the association between PM 10 and ischemic stroke attack. OR indicates odds ratio; CI, confidence interval. 
doi:1 0.1 371 /journal.pone.0095682.g007 



standardized measurement should be established. (3) In our 
research, meteorological condition was not included as a covariate, 
because only some of retrieved studies provided mean or median 
temperature data. Besides, the association between PM 10 and 
stroke attack was different from season to season [16]. Hence, it's 
not suitable to regard the mean or median temperature data as 
covariate in our research. Nevertheless, meteorological condition 
was an important covariate between studies. What's more, age is 
another important covariate to be considered, but current studies 
didn't provide sufficient data, they just provide the age range of 
study objects rather than the mean age. Future researches should 
provide more accurate data of age, temperature and pay more 
attention to the influence of age and temperature to stroke attack. 
(4) We only detect study design and area as important influence 



factors among PM and stroke by Meta-regression and subgroup 
analysis, however, there might exist other covariates that we didn't 
focus on, such as age, historical disease and meteorological 
condition. Besides, adding more studies may bring us more 
valuable covariates. (5) Lots of studies showed that different 
compositions of PM could cause diverse health outcomes. 
However, till 2012, there were no such articles, except Halonen's 
[23] study, which could provide detail data concentrating on the 
short-term effects of some composition of PM on stroke attack. In 
some studies [10-13,15,18-20], authors just mentioned traffic 
exhaust emission was the main origin of PM without providing the 
accurate data about the association between some composition of 
PM and stroke attack. While in other studies [7,8,16,21,22,24-26], 
authors even didn't mention the effect of composition of PM to 
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Table 2. Meta-analysis of subgroup studies. 







Covariate 


Level 


Number 
of articles 


Number 
of studies 


Model 


Pooled OR 

(95%CI) 


The effect of PM 2 . 5 
to stroke attack 


study design 


time-series study 


6 


22 


random effect model 


0.998 (0.993-1.002) 






case-crossover study 


3 


11 


fixed effect model 


1.029 (1.003-1.055) 


The effect of PM 10 
to stroke attack 


area 


Asia 


6 


17 


random effect model 


1.006 (1.000-1.011) 






Europe 


5 


39 


random effect model 


1.008 (1.002-1.014) 






North America 


2 


7 


fixed effect model 


1.001 (0.996-1.006) 






Oceania 


1 


5 


fixed effect model 


0.982 (0.974-0.990) 



PM 2 .5 indicates particular matter with aerodynamic diameter <2.5 um; PM 10 indicates particular matter with aerodynamic diameter <10 ujn; OR, odds ratio; CI, 

confidence interval. 
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stroke attack. As a result, we could only focus on the whole PM's 
effect on stroke attack. We urged future studies should pay more 
attention to the different composition of PM on health effect. 

Conclusions 

This Meta-regression and Meta-analysis raised significant issues 
that might help guide the future research in this area. PM2.5 and 
PM 1() exposure had different effect on different stroke subtypes. So 
it's worthwhile to study the effects of PM to ischemic stroke and 
hemorrhagic stroke, respectively. Standardizing of exposure 
measurement, bringing more covariates, exploring mechanisms 
and adding more studies in respective subtype are needed. 

Supporting Information 

File SI Table SI, sensitivity analysis of Meta-analysis in different 
particular matter to different stroke types. PM10 indicates 
particular matter with aerodynamic diameter <10 um; PM2.5, 
particular matter with aerodynamic diameter £2.5 (J,m; OR, odds 
ratio; CI, confidence interval. Table S2, meta-regression of 



different particular matters exposure to different stroke types. 
PM2.5 indicates particular matter with aerodynamic diameter £ 
2.5 |Xm; PM10, particular matter with aerodynamic diameter £ 
10 um; Coef, regression coefficient; Std. Err., standard error of 
logOR; CI, confidence interval. *Sex was not contained as a 
covariate for these studies were all about the whole population. 
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